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We consider a modification of the contact process incorporating 
higher-order reaction terms. The original contact process exhibits 
a non-equilibrium phase transition belonging to the universality 
class of directed percolation. The incorporated higher-order reaction 
terms lead to a non-trivial phase diagram. In particular, a line 
of continuous phase transitions is separated by a tricritical point 
from a line of discontinuous phase transitions. The corresponding 
tricritical scaling behavior is analyzed in detail, i.e., we determine 
the critical exponents, various universal scaling functions as well as 
universal amplitude combinations. 
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1 INTRODUCTION 

The concept of universality is well established for equilibrium critical phenomena 
where a unifying framework exists. Compared to the equilibrium situation less is 
known in case of non-equilibrium phase transitions. In particular a classification 
scheme of the rich and often surprising variety of non-equilibrium phase transitions 
is still lacking. Nevertheless, it is expected that analogous to equilibrium each non- 
equilibrium universality class is characterized by a certain symmetry. In this work 
we focus on stochastic processes which exhibit irreversible phase transitions into 
absorbing and fluctuation free states (see P] 01 El El for recent reviews) . In that 
case the universality determining symmetry is expressed in a corresponding path 
integral formulation [21 E] which is associated to the considered stochastic process. 
A well known example is the universality class of directed percolation (DP). Ac- 
cording to its robustness and ubiquity (including critical phenomena in physics, 
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biology, as well as catalytic chemical reactions) directed percolation is recognized 
as the paradigm of non-equilibrium phase transitions into absorbing states. All 
non-equilibrium critical systems belong to the DP universality class if the associ- 
ated absorbing phase transition is described by a single component order parameter 
and if the corresponding coarse grained system obeys the so-called rapidity rever- 
sal symmetry (at least asymptotically). Unfortunately, this symmetry is usally 
reflected on a coarse grained level only. Often, it is masked on a microscopic level, 
e.g. it is not reflected in the dynamic rules of certain lattice models. But the ro- 
bustness of the DP universality class is expressed by the conjecture of Janssen and 
Grassberger [E1EJ : short-range interacting systems, exhibiting a continuous phase 
transition into an absorbing state, belong to the DP universality class, if they are 
characterized by a one-component order parameter and no additional symmetries 
and no quenched disorder. This robustness of the DP universality class was recently 
demonstrated for five lattice models [S]. Despite the different interaction details 
(such as different lattice structures, different update schemes, different implementa- 
tion schemes of the conjugated field, infinite and finite numbers of absorbing states, 
models with and without multiple particle occupation, etc.) all considered models 
are characterized by the same universal scaling functions. Different universality 
classes than DP occur if the rapidity reversal is broken, e.g. by quenched disor- 
der 110[ I12[ I13[ I14j , or additional symmetries such as particle-hole symmetry 
(compact directed percolation |15[ I16j). particle conservation (Manna universal- 
ity class |17[ I18| I19j ) , or parity conservation (for example branching annihilating 
random walks with an even number of offsprings |2()1 121j ) . 

Another universality class of absorbing phase transitions which is directly re- 
lated to DP is tricritical directed percolation (TDP). Analogous to the c/> 6 -theory 
in equilibrium critical phenomena the process of TDP incorporates higher-order re- 
action terms than ordinary DP. Investigations of TDP in one-dimensional systems 
traces back to the seminal work of Grassberger [7j. Later, a field theoretical anal- 
ysis was performed by Ohtsuki and Keyes [HI EI] (see also 01231) • ^ n that work 
the authors discuss a one-component reaction-diffusion system that exhibits a tri- 
critical point. This tricritical point separates continuous DP-like transitions from 
first-order transitions. Using an e-expansion several critical exponents were esti- 
mated. Furthermore, the authors determine the upper critical dimension D c = 3. 
Compared to the established field theory of TDP less work was done numerically. 
Several modifications of one-dimensional DP-lattice models are known which yield 
multicritical behavior. But mostly a bicritical point is observed or the occurring 
tricriticality is analyzed within a mean field level only (see e.g. |25l 1261 l2"7j l. Sur- 
prisingly no systematic numerical investigations of tricritical-DP scaling behavior 
were performed so far. 

In this work we analyze the tricritical contact process (TCP). This model is a 
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modification of the well-known contact process |2H] which belongs to the universal- 
ity class of directed percolation. The modification incorporates higher-order parti- 
cle reactions, in particular a pair reaction scheme is added. It is worth mentioning 
that multicritical behavior does not necessarily appear by introducing additional 
higher-order reaction schemes. Multicritical behavior occurs if the resulting lower- 
order reactions vanish on a coarse grained scale (see [23] and references therein). 
In our case, the added pair reaction scheme leads to a non-trivial phase diagram 
containing first- and second-order transitions as well as a tricritical point. In D > 2 
the tricritical point separates a line of DP-like transitions from a line of first-order 
transitions. In one-dimension, preliminary results do indicate that the TCP does 
not exhibit a tricritical behavior. According to numerical simulations and to a 
mean field cluster analysis |29| the first-order line and the tricritical point collapse 
for D = 1 to the end point of the corresponding phase diagram (this agrees with 
the more general argument that first-order transitions can not occur in fluctuating 
one-dimensional systems due to the fact that the surface tension of a given domain 
does not depend on its size). Therefore, we focus on the higher-dimensional sys- 
tems in the following. First, we survey the mean field behavior of TDP and of 
the TCP which are valid for D > 3. Second, the two-dimensional TCP is numeri- 
cally examined. A scaling analysis of the DP-like transitions and of the tricritical 
phase transition itself is presented. In particular, the full crossover between both 
universality classes is recovered. The obtained values of universal quantities, such 
as the critical exponents, are compared to the results of the corresponding field 
theories 01221123]. 

Furthermore, the regime of first-order transitions is investigated. Discontinuous 
phase transitions are of their own interest (see for example |3U1 1261 1311 132j ) because 
most absorbing phase transitions are of second-order. As well known, discontinuous 
transitions are usually characterized by hysteresis cycles caused by the coexistence 
of two phases. A small but finite hysteresis can be observed in numerical simulations 
of the TCP. 

2 TRICRITICAL DIRECTED PERCOLATION 

The process of directed percolation might be represented by the Langevin equation 
(see e.g. 0) 

A" 1 d tPa = rp a - gp\ -cpl + T V V + h + rj. (1) 

The particle density on a mesoscopic (coarse grained) scale p a = p a (x,t) corre- 
sponds to the order parameter of the non-equilibrium phase transition. The control 
parameter of the transition r describes the distance from the critical point r = 0. 
A finite positive particle density occurs above the transition point (r > 0) whereas 
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the absorbing vacuum state (p a = 0) is approached below the transition point. The 
external field h is conjugated to the order parameter and is usually implemented 
as a spontaneous particle creation process (see e.g. Furthermore, r/ = rj(x,t) 

denotes the noise which accounts for fluctuations of the particle density. The noise 
r\ is a Gaussian random variable with zero mean and whose correlator is given by [H] 

( v(x, t)v{x'i 0) = ^~ lft PMi t) % - x') - *')• ( 2 ) 

Notice, that the noise ensures that the systems is trapped in the absorbing state 
Pa,(x,t)=0. Furthermore, higher-order terms such as /Q a (ai, t) 3 , /0 a (x, t) 4 , . . . or 
V 4 /0 a (x, t), V 6 p a (x, t), . . . are usually neglected because they are irrelevant under 
renormalization group transformations as long as g > 0. Negative values of g give 
rise to a first-order phase transition whereas g = is associated with a tricritical 
point j2HHsIl- I n the latter cases the cubic term p^ cannot be neglected. 

In the following we present a simple but instructive mean field treatment. Ne- 
glecting the noise term as well as spatial variations of the order parameter the 
steady state behavior (<9 t/ o a = 0) at zero field is given by 




The first solution corresponds to the absorbing phase. According to a linear sta- 
bility analysis it is stable for r < and unstable for (r > 0). The solution with 
the — sign yields unphysical and unstable results. The +sign solution describes 
the order parameter as a function of the control parameter r and of the additional 
scaling field g. For g < 0, this solution is stable if t > —g 2 /Ac, otherwise it is 
unstable. Assuming that the system is in the active phase (/9 a > 0) the order pa- 
rameter jumps at the borderline (r = —g 2 /4c) from p a = \g\/2c to zero. Thus the 
absorbing and the active phase coexist between the spinodal r = — g 2 /Ac and the 
line r = for g < 0. 

The tricritical behavior is obtained for g = 0. The order parameter behaves for 
r > as 

Pa = (r/ct (4) 

with the tricritical order parameter exponent (3 t = 1/2. For g > 0, the active phase 
is stable for r > and unstable otherwise. Close to the transition points (r = 0) 
the order parameter exhibits a directed percolation like behavior 

- +0(t 2 ) (5) 
9 J 

with the exponent /? DP = 1. The amplitude of this power-law diverges for g — > 0, 
signaling the crossover to the tricritical behavior. 
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Figure 1: The mean field phase diagrams of tricritical directed percolation (TDP). 
The left figure sketches the phase diagram as a function of the coarse grained 
variables r and g [see Eq. (fT])]. The right figure shows the corresponding phase 
diagram as a function of the microscopic variables p and q, as defined by a reaction 
scheme of the tricritical contact process (see text). The thick lines correspond 
to second-order phase transitions (Il.order) with /3 DP = 1. The bold circles mark 
the transition point of TDP with j3 t = 1/2. The shadowed areas indicate the 
coexistence of the absorbing and the active phase. The thin dashed lines illustrate 
the crossover to the tricritical behavior. In the right figure, two absorbing phases 
appear for q = 1 (long dashed lines). These phases are the fully occupied lattice 
(stable for q > 0) and the empty lattice (stable for q < 1/2). 



Sufficiently away from the critical line r = we observe again the tricritical 
behavior 



Pa 



2c 



2c J 



for both g > and g < 0. The crossover from the tricritical behavior takes place 
at 

r = 0(g 2 /4c) . (7) 

The complete phase diagram is sketched in Fig.^ 

Within the active phase for g > the order parameter obeys the scaling form 



P, = \- pt r tDP (Ar, g\t h = 0) 



(8) 



for all positive A and with the crossover exponent <j) = 1/2- In particular, A = g 1 ^ 
leads to the scaling form 



(9) 
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with 

{4 P if x <C 1 
(10) 
x Pt if x > 1 . 

In that way, the scaling plot p a ./g^ t ^ vs. rg^ 1 ^ reflects the crossover from DP to 
TDP. 

Now we consider the scaling behavior at non-zero conjugated field h. The 
tricritical order parameter behavior at criticality (r = and g = 0) is given by 

Pa = (h/c)M« (11) 

with the exponent <r t = 3/2. At the tricritical point the equation of state obeys 
the scaling form 

p a = A" A f tDP (Ar, g = 0, h\^) . (12) 
Instead of the above equation the so-called Widom-Griffiths scaling form 

h = \-^ h tDP (\r, p a \&,g = 0) (13) 

is often used to describe the scaling behavior of the equation of state. Despite 
metric factors the tricritical scaling form is given by 

ht OP (x,y,g = 0) = -xy + y 3 . (14) 

The corresponding scaling form of ordinary DP h DP (x,y) = —xy + y 2 reflects the 
different universality class. 

For the sake of completeness we present the dynamical order parameter behavior 
at tricriticality 

p a (r = 0, g = 0, h = 0, t) = (p~l + 2ct ) 

(2ct)~ at (15) 

t— >oo 

with o t = 1/2. Furthermore, the steady state value of the order parameter [Eq. 
is approached from below as 

p a (r, g = 0, h = 0, t) = ^ [l - Coe-^l + 0(e~ 2t ^ )] . (16) 

Here, the constant c contains the initial conditions and £y denotes the temporal 
correlation length 

e„ = (2t)-»U (17) 

with z^t = 1- Incorporating spatial variations of the order parameter the spatial 
correlation length £ ± can be derived via an Ornstein-Zernicke approach. The re- 
sulting correlation length exponent z/ ± t = 1/2 leads to the dynamical exponent 
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z t = V\\ :t /ux,t = 2. Eventually we just mention that the tricritical exponent of 
the order parameter fluctuations is given by 7^ t = 1/2 [221 ■ This value reflects a 
qualitative difference between the mean field scaling behavior of DP and TDP. In 
the latter case the fluctuations diverge at the transition point whereas they remain 
finite (jump) in case of DP |341 I35| . Another difference between both universal- 
ity classes concerns the value of the upper critical dimension, namely D c = 4 for 
ordinary DP [SHI EI! and D c = 3 for TDP (211231. 



3 MODEL AND SIMULATIONS 

We consider a modification of the contact process (CP). The CP was introduced 
by Harris in order to model the spreading of epidemics |28j . It is a continuous-time 
Markov process that is usually defined on a D-dimensional simple cubic lattice. 
A lattice site may be empty (n = 0) or occupied (n = 1) by a particle and the 
dynamics is characterized by spontaneously occurring processes, taking place with 
certain transition rates. In numerical simulations the asynchronous update is re- 
alized by a random sequential update scheme: A particle on a randomly selected 
lattice site i is annihilated with rate one, whereas particle creation takes places on 
an empty neighboring site with rate XN/2D, i.e., 

n, = l -j* n l = 0, (18) 

Tli = — ► Hi = 1 , (19) 
XN/2D ' V 7 

where N denotes the number of occupied neighbors of n 4 . Note that the rates 
are defined as transition probabilities per time unit, i.e., they may be larger than 
one. Thus, rescaling the time will change the transition rates. In simulations 
a discrete time formulation of the contact process is performed. In that case a 
particle creation takes place at a randomly chosen neighbor site with probability 
p = A/(l+A) whereas particle annihilation occurs with probability 1—p = 1/(1+A). 

In the language of reaction-diffusion models the contact process is described by 
the reduced reaction scheme 



A — > , A — > 2A (20) 

which is controlled by the parameter p and where the quantity A represents a 
particle. According to Ohtsuki and Keyes |22U23j higher-order reactions may lead 
to a tricritical behavior as well as to a first-order behavior. Here, we additionally 
take into consideration the pair reaction 



2A — ► 3A 



(21) 
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This reaction is controlled by a parameter q. Updating a given occupied lattice 
site i, we first perform the pair reaction scheme with probability q (otherwise a usual 
CP-update step is performed). If a randomly selected neighbor of i is occupied we 
add a third particle at an empty (also randomly selected) neighbor of the pair. If 
the lattice site i is isolated a usual update procedure of the CP is performed. The 
detailed reaction scheme is listed in Tabled 

It is essential for the following that within our implementation the annihilation 
processes 

A — > , AA — > OA , AA — ► AO (22) 

take place with probabilities proportional to (1 — q). Thus these reactions are 
suppressed for q ~ 1 and are eventually forbidden for q = 1. In that case particle 
annihilation occur via the reactions 

AO — ► 00 , 0,4 — > 00 . (23) 

This implies that annihilation processes are restricted to the domain boundaries 
and do not occur insight of a domain. This behavior is reminiscent of compact di- 
rected percolation ^HHEIEII where the density of particles exhibits a discontinuous 
behavior at the critical point. 

The above scenario is also reflected by a corresponding mean field analysis. 
Within the simplest mean field approach (single site approximation), the order 
parameter is given by 

dtp, = r Piq p a - g Ptq pl - c p>q p\ (24) 

with 

T P , q = 2p-l, g p>q = p - (2 - p)q , c p , 9 = q . (25) 

The tricritical point (p t = 1/2, q t = 1/3) is determined by the conditions r = 
and g = 0. The line of second-order phase transitions (p c = 1/2 and q < 1/3) is 
obtained from r = and g > 0. Within the first-order regime (g < 0) the absorbing 

Table 1: The reaction schemes and the corresponding probabilities of the con- 
sidered tricritical contact process. For q = the reaction scheme of the original 
contact process is recovered (see text). 



reaction 


probability 


A — ► 
0.4 — ► 00 
OA — >AA 
AA — ► 0^ 
0AA — ► AAA 


(l-g)(l-p) 
g(i-p) 

(1 - Q)P 
(1-?)(1-P) 

q 
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phase is stable for p < 1/2. The borderline of stability of the active phase (p a > 0) 
is determined by r = —g 2 /4c. The corresponding mean field phase diagram is 
sketched in Fig.[I] 

For q = 1 the corresponding mean field differential equation is given by 

cVa = (2p - l)p a (l - p a ) + pl(l - p a ) . (26) 

Obviously, the steady state solutions are the empty lattice (p a = stable for 
p < 1/2) and the fully occupied lattice (p a = 1 stable for p > 0). Note that both 
phases coexist for p < 1/2. It is worth to compare the TCP for q = 1 to the process 
of compact directed percolation which is described by the equation (see e.g. (Sj) 

d t p a = (2p-l)p a (l-p a ). (27) 

Here, the fully occupied lattice is stable for p > 1/2 whereas the empty lattice is 
stable for p < 1/2. Thus the process of compact directed percolation displays no 
phase coexistence in contrast to the tricritical contact process for q = 1. Further- 
more, Eq. (|27|) clearly express the particle-hole symmetry which is the characteristic 
symmetry of the universality class of compact directed percolation. In case of the 
TCP for q = 1, the particle- hole symmetry is broken by the pair reaction processes 
contributing p 2 (l — p a ) to the equation of motion. 

Analyzing numerically or experimentally the scaling behavior of tricritical sys- 
tems it is crucial to determine the value of the tricritical point with high accuracy. 
Thus a sensitive criterion is required to distinguish a first-order transition from a 
second-order transition. At first glance, one is tempted to make use of the order 
parameter jump at the first-order transition. At the tricritical point the order pa- 
rameter changes its behavior from a discontinuous jump to a continuous power-law. 
But first, this behavior is affected close to the tricritical point by crossover effects 
(see Fig.^). Second, it is notoriously difficult to distinguish a continuous phase 
transition with a small but finite value of the exponent (3 t from a discontinuous 
jump. Especially this situation occurs in case of two-dimensional TDP. An alter- 
native way is to investigate instead of the order parameter the order parameter 
fluctuations which diverge at the tricritical point but remain finite within the first- 
order regime. But the fluctuation measurements suffer by crossover effects in a 
similar way as the order parameter, i.e., a corresponding analysis yields therefore 
no significant improvement. 

In order to circumvent these problems we apply a method of analyzing that is 
based on the scaling behavior of the order parameter close to the tricritical point 
within the second-order transition regime. In particular, the scaling form of the 
order parameter is used to recover the complete crossover from ordinary DP to 
the TDP. We assume that the order parameter as well as the order parameter 
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fluctuations obey the scaling forms 

Pa ~ A~ ft f tDP (Ar, gX*,h = 0) (28) 

Ap a ~ A 7 * d tDP (Ar, g\t h = 0) (29) 

with the so far unknown tricritical exponents (3 t , j' t , and <j>. Asymptotically, the 
scaling functions have to fulfill the power-laws 

x Pdp jf x <C 1 
f tDP (x,l,0) ~ { (30) 
x Px if x > 1 , 

x _7 dp if j;<1 
d tDP (x,l,0) ~ { (31) 
x 7t if 

The values of the exponents /3 DP and 7^ p are known with sufficient accuracy (see 
TableEl and references therein). A serious problem is caused by the fact that 
the scaling forms contain the coarse grained variables r and g. In general, these 
variables depend on the (microscopic) parameters p and q in an unknown way. But 
there exists a certain window of scaling where the coarse grained variables can be 
replaced by the model parameters p and q. This will allow the determination of 
the tricritical point as well as of the exponents j3 t and 7'. 

In the following the system is investigated within the active phase close to the 
line of second-order phase transitions, i.e., g(p,q) > and r(p,q) close to zero. 
Performing simulations, the order parameter is measured as a function of p for 
fixed q. Sufficiently close to the transition line p c (q) the coarse grained variable is 
approximated 

r(p,q) = T(p c (q) +Sp,q) 

dr 



=0 



5p (32) 

Pc(q) 



with bp = p — p c (q). Since the derivative is taken at the transition line p c (q) its 
value depends on the parameter q. But close to the tricritical point we yield up to 
higher-orders 

r( P ,q) - I 
with Sq = q t — q. 

A similar approximation can be obtained for g(p,q). In the vicinity of the 
critical line we use the approximation 



Sp + 0{5p5q) (33) 

Pc(gt) 



g(p,q) = g(pc(q) + Sp,q) 

~ g(Pc(q),q) + ^ 



5p . (34) 

Pc(q) 
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Figure 2: The order parameter p a of the two-dimensional contact process (CP) as 
a function of the control parameter p. The data are obtained from simulations on a 
simple cubic lattices of linear size L for various field values (from 10 -6 to h = 1CP 5 ). 
For non-zero field p a exhibits an analytical behavior (lines). For zero field (symbols) 
the order parameter vanishes at the transition point p = 0.62246 |39j . 

Sufficiently close to the tricritical point we obtain 

g(p,q) « Vg-t 5q 
pt,qt 



+ dg 
dp 



5p + 0(5 P 5q), (35) 



pt,gt 



with Vg = (d p g, d q g) and where the tangential vector along the phase boundary is 
denoted by t = (dp c (q)/dq, 1). Assuming that Vg ■ t and dg/dp are of the same 
order at the tricritical point, the coarse grained variable g(p, q) can be replaced in 
the scaling functions Eqs. (JHUl ITTT^ by the reduced model parameter 5q if 

5p < 5q . (36) 

This condition Eq. (|36|) is fulfilled if the simulations are performed in a way that 
the distance to the phase boundary is smaller than the distance to the tricritical 
point. Thus, the order parameter and its fluctuations obey for 5q > the scaling 
forms 

Pa ~ A~ ft r q (\Sp, 6q\*, h = 0) , (37) 

Ap a ~ A 7 * d q (\5p, <5gA^, h = 0) (38) 

if the critical point p c (q) is approached along q = const paths (as indicated by 
the index q). The above scaling forms are not valid if _L t or if the derivative 
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drjdp vanishes at the tricritical point. In the latter case higher-orders (0(Sp 2 )) 
are required to describe the scaling behavior. 

In the next section we describe the analysis of the simulation data. In partic- 
ular, we have measured the averaged density of active sites p^[p,q,h) = {L~ D N a ), 
i.e., the order parameter as a function of p and q and of the conjugated field h. 
Numerically obtained order parameter curves for q = and for different values of 
the conjugated field are plotted in Fig.|2l The conjugated field is implemented via 
a generation of particles (0 — > A with probability h). It results in a rounding of 
the zero field curves and the order parameter behaves smoothly as a function of the 
control parameter for finite field values. For h — > we recover the non-analytical 
order parameter behavior. Additionally to the order parameter, we investigate the 
order parameter fluctuations Ap a (p, g,h) = L D ((pl) — (p a ) 2 ) and the susceptibil- 
ity xiPi <7) h) = dp a /dh. The susceptibility is obtained by performing the numerical 
derivative of the order parameter p a with respect to the conjugated field h. 

4 SECOND-ORDER PHASE TRANSITION: 
DIRECTED PERCOLATION BEHAVIOR 

According to the above analysis we have simulated the tricritical contact process 
(TCP) and have measured p a as a function of p, keeping q fixed. Order parameter 
curves for various g-value are shown in Fig.El According to Eq. (|H7|) these different 

1.0 
0.8 
0.6 
0.4 
0.2 



"0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
P 

Figure 3: The order parameter p a of the tricritical contact process (TCP) as a 
function of the control parameter p for various values of q G {0, 0.1, 0.2, 0.3, . . . , 0.9} 
(from right to left). The dashed line corresponds to the pure contact process 
(q = 0). The data are obtained from simulations on simple cubic lattices of linear 
size L = 64,128, ... ,512. 
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Figure 4: The crossover scaling function of the order parameter at zero field. The 
dashed lines correspond to the asymptotic behaviors, i.e., to the ordinary DP and 
to the tricritical DP behavior [see Eq. . 

curves collapse onto a single curve if the rescaled order parameter p a 5q~^^ is 
plotted as a function of the rescaled control parameter 5p5q~ l l^. Therefore, we 
vary the parameters /3 t , </> as well as q t until a data collapse is obtained. Convincing 
results are obtained for (3 t = 0.14 ± 0.02, 4> = 0.55 ± 0.03 and q t = 0.9055 ± 
0.0020. A corresponding scaling plot is shown in Fig.0] The data are obtained 
from simulations for 16 different values of q ranging from q = 0.55 up to q = 0.904. 
Typical distances to the critical line are of the order of 0(5p) = 10 -4 . On the 
other hand the minimal distance Sq to the tricritical point along the g-axis is larger 
than 0.0015. In that way, the condition Eq. (|36|) is fulfilled, justifying the use of 
the scaling forms. 

Since the entire crossover region covered several decades it could be difficult to 
observe small but systematic differences between the scaling functions of both mod- 
els. It its therefore instructive to examine the crossover via the so-called effective 
exponent 

A* = t^— lnf,(x,l,0). (39) 

The corresponding data are shown in Fig.[5J The excellent data collapse of /3 e fj over 
more than 6 decades reflects the accuracy of the determination of the tricritical 
point. 

A similar crossover scaling analysis can be performed for the order parameter 
fluctuations. In that case, we use the above determined values of /3 t , q t and vary the 
value of the fluctuation exponent j' t until a data collapse of the different g-curves 
occurs. Figure^] shows the corresponding data collapse as well as the effective 
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Figure 5: The effective exponent /9 eff of the order parameter. Both asymptotic 
scaling regimes (f3 t and /3 DP ) as well as the crossover regime are clearly recovered. 

exponent 

= ^ln^(ar,l,0). (40) 

Although the data of the effective exponent are suffering from statistical fluctu- 
ations both asymptotic scaling regimes as well as the crossover regime can be 
identified. Worth mentioning, the crossover scaling function exhibits a so-called 
non-monotonic crossover |40| . Unfortunately, the statistical scattering of the ef- 
fective exponent data masks the non-monotonic behavior from j' t = 0.93 ± 0.06 
to 7^ P . 

In summary, the crossover scaling analysis allows the determination of the 
tricritical value q t with high accuracy. Furthermore, the effective exponents f3 eS 
and 7' ff reflect the full crossover from tricritical DP to ordinary DP which spans 
more than 6 decades. 

5 SECOND-ORDER PHASE TRANSITION: 
TRICRITICAL BEHAVIOR 

In the previous section we have determined the tricritical value q t as well as the 
critical line q(p c ) of the directed percolation like phase transitions. This yield a 
first glance of the phase diagram which is presented in Fig.0 We now investigate 
the tricritical scaling behavior in detail. Therefore, we consider first the order 
parameter at zero field in the vicinity of the tricritical point. Second a conjugated 
field is applied which allows the determination of various universal quantities such 
as tricritical exponents, scaling functions as well as amplitude combinations. 
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Figure 6: The crossover scaling function of the order parameter fluctuations at 
zero field. The dashed lines correspond to the asymptotic behaviors, i.e., to the 
ordinary DP and to the tricritical DP behavior [see Eq. (|31|) ]. The inset displays 
the corresponding effective exponent j' eS . 



In the following analysis, the tricritical point is approached in three different 
ways and the order parameter is determined at zero field (see Fig.EJ). First, p a 
is examined within the active phase as a function of p at the tricritical value 
q = 0.9055. This yields the estimate 

Pt = Pc (q t = 0.9055) = 0.29931 ± 0.00003 . (41) 

Second, the order parameter is measured for p = 0.29931 as a function of q. From 
this analysis we obtain 

q c (p t = 0.29931) = 0.90552 ± 0.00005 (42) 

which agrees with the above determined value q t = 0.9055 ±0.0020. To check these 
estimates the critical point is approached along a third path which is (more or 
less) perpendicular to the phase boundary. As well known the leading order of the 
scaling behavior, i.e., the critical exponents, do not depend on the way the critical 
point is approached. But the prefactors of the corresponding power-laws and the 
corrections to the leading scaling order are affected by the different directions. Thus 
the scaling behavior of the order parameter obeys asymptotically 

p a ~ (dpSpf 1 for q = q t (43) 
p a ~ [a q 5q) Pt for p = p t (44) 
p a ~ (a ± 5p)^ 1 along the perpendicular path. (45) 
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Figure 7: The phase diagram of the two-dimensional modified contact process. The 
solid line marks continuous phase transitions which belong to the universality class 
of directed percolation. The dashed line corresponds to first-order phase transitions 
and the bold circle indicates the tricritical point. The small circles in the right 
figure show where the transitions points are determined numerically. Furthermore, 
the tricritical point is approached in the simulations along three different ways 
illustrated by the three arrows. 

Plotting the corresponding data accordingly (e.g. p a as a function of (a p 5p)^*) 
the three different curves collapse asymptotically to the leading tricritical behavior 
if the so-called non-universal metric factors (Zp, & _j_ £11*6 chosen appropriately. 
The leading tricritical behavior corresponds in that analysis to a straight line with 
slope one. As can be clearly seen in Fig. |H] all three curves approach the tricritical 
scaling behavior asymptotically, confirming the accuracy of the determination of 
the tricritical point and of the order parameter exponent j3 t . Furthermore, the 
corrections to scaling, i.e., the deviations to the asymptotical power-law depend 
strongly on the way the tricritical point is approached. Surprisingly, the smallest 
corrections occurs along the path p = const and not as usually expected along the 
way which is perpendicular to the phase boundary. 

So far we have considered the order parameter at zero field in the vicinity of 
the tricritical point. Applying an external field which is conjugated to the order 
parameter, it is possible to investigate the scaling behavior of the tricritical equation 
of state. Again, the conjugated field is implemented via a generation of particles 
(0 — > A with probability h). Sufficiently close to the tricritical point (g = 0) the 
order parameter obeys the scaling form 

P,{r P , q ,g P , q ,h) ~ A _/3t r(Ar p .„ gX*, hX^ ) 

3=0 

= X-^ f(Xr p , q ,0,hX^). (46) 

Here, r Pi9 describes the distance to the tricritical point within the p — q— plane. 
All non-universal system dependent features (such as the lattice structure or the 
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Figure 8: The order parameter behavior of the tricritical contact process (TCP). 
The tricritical point is approached along three different ways indicated by the three 
different symbols. All three curves tend asymptotically to the function f(x) = x 
(the dashed line corresponds to the pure tricritical power-law behavior) if p a is 
plotted as a function of (a p 5p)^ , (a^q)^, and (a±6p)^ t , respectively. But as can 
be seen the corrections to scaling depend on the scaling direction. 

way the tricritical point is approached, etc.) can be absorbed in two metric factors 
dpath and a h . Once the non-universal metric factors are chosen in a specific way, 
the scaling function R (in contrast to f) is the same for all systems belonging to 
the universality class of directed percolation. Thus the universal scaling behavior 
is described by the ansatz 

P,(r P , q ,0,h) ~ i?(Aa path r path ,0,a h /iA fft ) (47) 

yielding for a h h\ at = 1 

p a (r Pig , 0, h) ~ (a h h)-^ J R(a path r path (a h / l ) 1 / <Tt , 0, 1) . (48) 

Throughout this work we use the conditions ii(l,0,0) = i?(0,0, 1) = 1 to specify 
the metric factors which can be obtained from the amplitudes of the corresponding 
power-laws [see e.g. Eqs. (|43 |l -(|45 |) ]. 

In order to determine the universal scaling form of the equation of state we 
have measured the field dependence of the order parameter. Again the tricritical 
point is crossed along three different paths in the p — q— plane. Along each path, 
p a is determined as a function of the distance to the tricritical point for at least 
four different field values. The results are plotted for the perpendicular path in 
Fig.|SJ According to the scaling form Eq. (|48[) the rescaled order parameter is 
plotted as a function of the rescaled distances to the tricritical point. In our 
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Figure 9: The order parameter p a of the two-dimensional tricritical contact process 
(TCP) close to the tricritical point. The tricritical point is crossed perpendicular to 
the phase boundary (see Fig^). The order parameter is plotted as a function of the 
control parameter p. The zero field data (circles) are obtained from simulations 
on a simple cubic lattices of linear size L = 64, 128, 256, 512. The solid curves 
correspond to the non-zero field behavior (from h = 2 10" 3 toh = 3 1CT 5 ). 



analysis we have varied the field exponent <r t until a data collapse is obtained. 
Convincing results are obtained for cr t = 1.20 ± 0.05 and are shown in Fig.1101 
Since data of different scaling directions are considered the presented data collapse 
is an impressive demonstration of the universality class of TDP. Furthermore, the 
tricritical universal scaling function differs significantly from the corresponding 
scaling function of ordinary DP, reflecting the different universality classes. 

Additionally to the order parameter behavior we have investigated the order 
parameter fluctuations A/j a as well as the order parameter susceptibility x- The 
scaling behavior of both quantities is described by [H] 

aAAp a (r M ,0,/i) ~ A 7 ' Z)(Aa path r path ,0,a h /iA fTt ), (49) 

a x x(T P , g ,0, h) ~ A 7t X(Aa path r path ,0,a h /iA CTt ). (50) 

Setting D(0, 0, 1) = 1 the non-universal metric factor a A is specified. Taking into 
account that the susceptibility is defined as the derivative of the order parameter 
with respect to the conjugated field we find a x = 

X(x,0,y) = d y R(x,0,y), (51) 
X(0,0,1) = A, (52) 
7t = <7 t - A ■ (53) 
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Figure 10: The universal scaling function R(x, 0, 1) of the universality class of TDP. 
For all three scaling directions (only the data within the scaling regime are plotted) 
at least four different curves are plotted corresponding to four different field values. 
The circle marks the condition i?(0, 0, 1) = 1. The dashed line corresponds to the 
universal scaling function of ordinary DP (taken from 8 ). 



The latter equation corresponds to the Widom scaling law in equilibrium. Further- 
more, Eq. (|52|) offers a useful consistency check of the numerical analysis. 

The universal scaling functions of the order parameter fluctuations and the or- 
der parameter susceptibility are shown in Fig.^J an d Fig. 1121 respectively. The 
susceptibility is obtained by performing the numerical derivative of the order pa- 
rameter with respect to the conjugated field. Both scaling functions exhibit a 
maximum signaling the divergence of A/9 a and x at the tricritical point. The sus- 
ceptibility data fulfill Eq. (|52|). reflecting the accuracy of the performed analysis. 
In both cases, the obtained universal scaling functions differ significantly from the 
corresponding scaling functions of ordinary DP. 

Additionally to the critical exponents and universal scaling functions it is useful 
to investigate universal amplitude combinations. An often considered amplitude 
combination is related to the susceptibility behavior below and above the transition. 
For e.g. q = const the susceptibility diverges as 

X (5p> 0,g = 0,h = 0) ~ a Xt+ 6p-^, (54) 

X (5p <0,g = 0,h = 0) ~ a x ,_ (-5p)~^ , (55) 

if the critical point is approached from above and below, respectively. Using 
Eq. (|50|) the susceptibility ratio 

X(5p> 0,0, h) X( a p 5p X,0,a h hX at ) 

X (5p< 0,0, h) ~ X(-a p 5p\,0,a h hX^) ... 
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Figure 11: The universal scaling function D(x, 0, 1) of the universality class of 
tricritical directed percolation. For all three scaling directions (only the data within 
the scaling regime are plotted) the scaling plot contains at least four different 
curves corresponding to four different field values. The circle marks the condition 
£(0,0,1) = 1 

= (56 ) 
X(-l,0,z) 

is clearly a universal quantity for all values of the scaling argument x = a h h\a p 8p\~ at . 
In particular it equals the ratio a Xi+ /a Xi _ for vanishing field 

«x, + = X(+1,0,0) 
a x ,_ X(-1,0,0) ' 

For example, the mean field value of the universal ratio equals 1/2. 

The ratio X(l, 0, x)/X(— 1, 0, x) is shown in the inset of Fig. 1121 The extrapo- 
lation to the tricritical point (x — > 0) yields the value of the universal amplitude 
ratio X(l,0,0)/X(-1,0,0) = 0.35 ± 0.05. This value differs significantly from the 
corresponding values of ordinary DP for D = 2 and D = 3 (see Fig.^J as well 
as Table[5J). Worth mentioning, the universal susceptibility ratio exhibits a non- 
monotonic behavior as a function of the scaling argument x. This non-monotonic 
behavior is a characteristic feature of TDP and does not occur in ordinary DP (see 
e.g. Fig. 30 of 0). 

So far simulation data are taken into account where the correlation length is 
small compared to the system size L. Thus, the data presented above do not 
suffer from finite-size effects, such as rounding and shifting of the anomalies. A 
typical feature of finite-size effects in equilibrium is that a given system may pass 
within the simulations from one phase to the other. This behavior is caused by 
critical fluctuations which increase if one approaches the transition point. In case 
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Figure 12: The universal scaling function X(x,0, 1) of the universality class of 
tricritical directed percolation. The circle marks the condition X(0, 1,0) = (3 t /&t 
and reflects the accuracy of the performed analysis. The inset displays the uni- 
versal ratio X(l, 0, x)/X(— 1, 0, x). The extrapolation (x — > 0) yields the value of 
the universal amplitude ratio X(l, 0, 0)/X(-l, 0, 0) = 0.35 ± 0.04. The obtained 
value differs significantly from the corresponding values (dashed lines) of ordinary 
directed percolation 0.25 for D = 2 and 0.65 for D = 3 [HE]; respectively. 



of absorbing phase transitions the scenario is different. Approaching the transition 
point, the (spatial) correlation length £ ± increases. As soon as £ ± is of the order 
of L the system may pass to the absorbing state where it is trapped forever. As 
pointed out in |3J an appropriate way to handle that problem is to incorporate a 
conjugated field. Due to the conjugated field the system can not be trapped within 
the absorbing state and steady state quantities are available for all values of the 
control parameter. 

As usual the system size L enters the scaling forms as an additional scaling 
field, e.g. 

p a - A _/3t £(Aa path r path , 0, a h h\°\ a L L\- u ^) (58) 

with the tricritical exponent u ± t of the spatial correlation length. Using the above 
scaling form it is possible to determine the correlation length exponent from data of 
different system sizes. As well known from equilibrium, ratios of order parameter 
moments (p^) are more suited to estimate the correlation length exponent. For 
example, the well-known Binder cumulant Q = 1 — (pf)/3(p a ) was successfully 
investigated in numerous works (see e.g. [12]) dealing with equilibrium as well as 
non-equilibrium phase transitions. Furthermore, the value of the Binder cumulant 
at the transition point is a universal quantity. 

Unfortunately, the Binder cumulant diverges at the critical point of absorbing 
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Figure 13: Scaling plot of the ratio U(h, L) close to the tricritical point (p = 
0.29931,(7 = 0.90552). For vanishing field the ratio tends to the universal value U = 
0.84 ± 0.04 which characterizes the universality class. The dashed lines correspond 
to the values of U(h — > 0) for two-dimensional and three-dimensional directed 
percolation. The inset displays the ratio slightly away from the critical point (p = 
0.29934,5 = 0.90557). The systematic deviations from a data collapse for small 
values of the scaling argument reflect how sensible U depends on the determination 
of the critical point. 




phase transitions |41[ 143 j. This behavior is caused by the vanishing steady state 
fluctuations in the absorbing phase and reflects the different nature of the zero- 
order parameter phase in equilibrium and in absorbing phase transitions. A ratio 
that remains finite at criticality is given by [13] 

(pD(pI) - (p,)(p 2 ) 2 

(Pa)(pf) " (P,)(P 2 J< 

This ratio is as useful for absorbing phase transitions as the Binder cumulant Q 
is for equilibrium, i.e., its value at criticality characterizes the universality class. 
Here, we investigated the ratio U close to the tricritical point (r = and g = 0). 
Its scaling behavior obeys 



U = /„ W„4X r, X / n2 \2 ' ( 59 ) 



Pt,9t 



U = U(\r p , q ,g p , q \' i> ,h\ a \L\- u 

= U(0,0,hX^,LX- v ^)\ LX _ VXit=i 

= U (0,0, hL^I^). (60) 

The ratio U is shown in Fig.^Jfor p = 0.29931 and q = 0.90552. Convincing data 
collapses are obtained for u±_ t = 0.59±0.08. As can be seen the ratio tends to a well 
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defined value for h — > independent of the system size L. The obtained value U = 
0.84 ± 0.04 differs significantly from the corresponding values of two-dimensional 
and three-dimensional ordinary DP (see Table|2J). Note that the determination 
of U(h — > 0) does not depend on the critical exponents a t and u ±>t . But it is 
very sensitive to the determination of the critical point. Performing simulations 
slightly away from the critical point the ratio U(h — > 0) displays a clear system 
size dependence (see inset of Fig. lT3*|) . 

Thus we have determined the steady state scaling behavior of tricritical directed 
percolation. The obtained values of the critical exponents as well as of the universal 
amplitude ratios are listed in Table[2j The accuracy of the estimated exponents 
can be checked with the scaling law 

7t ' = u Xit D - 2/? ±it . (61) 

The determined two-dimensional values 7' = 0.93 ± 0.06, u ±tt = 0.59 ± 0.08, and 
(3 t = 0.14 db 0.02 fulfill the above scaling law within the error bars. Furthermore, 
the exponent of the spatial correlation function rj± :t is related to the correlation 
length exponent and to the fluctuation exponent via the Fisher scaling law 

(2 - ru.tKt = l[ , (62) 

leading to r] ± ^ t = 0.42±0.24. Another quantity of interest is the fractal dimension D f 
of growing clusters at criticality. The fractal dimension is given by |44| 

D f = D — (63) 

yielding D { = 1.76 ± 0.05. This value is larger than the corresponding values 
of ordinary directed percolation (D { D=2 ~ 1-20 and D f D=3 w 1.60). A detailed 
investigation of the fractal behavior of critical clusters is desirable. For example, a 
determination of the lacunarity along the phase boundary would present a deeper 
understanding of the cluster propagation (see e.g. [Tl 

It is worth comparing our numerical results to those of corresponding field 
theoretical analyses. Within a two-loop approach the critical exponents are given 
in linear order of e = £>„ - D bv [221 1231 IB] 

A = - - e0.4580... , (64) 

P't = l + 0(e 2 ) (65) 
z t = 2 + e 0.0086... , (66) 

u ±it = - + e 0.0075... , (67) 
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Table 2: The critical exponents and various universal amplitude combinations of 
tricritical (TDP) and ordinary directed percolation (DP) for various dimensions D. 
The crossover exponent from tricritical to ordinary DP is given by 4> D=2 = 0.55 ± 
0.03 and (f) MF = 1/2. 
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0.14 ±0.02 
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X(-1,0) 


0.35 ±0.05 
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0.25 ±0.01 


0.65 ±0.03 
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0.84 ±0.04 




0.704 ± 0.013 


0.61 ±0.02 


1/2 



7t = l+0(e 2 ) (68) 
7t ' = \ + e 0.4386... , (69) 

0=2 - e0 - 0121 -- - • ( 70 ) 

Here, (3' denotes the critical exponent of the survival probability. According to the 
e-expansion one expects that the two-dimensional (e = 1) values of the exponents 
Vx,ti 7t; an d differ only slightly from their mean field values. Whereas strong 
deviations are predicted by the e-expansion for (3 t and r y[. This behavior is con- 
firmed by our numerical results. But one has to mention that the field theorically 
estimated exponents differ significantly from the numerical values. For example, 
the e-expansion yields for the order parameter exponent f3 t = 0.042 .... This value 
differs by 70% from the numerical value. Furthermore, the RG results predict 
4> < 1/2 whereas the simulations clearly show 4> > 1/2. Thus an e-expansion of 
higher-orders than 0(e) is desirable to describe the scaling behavior by a field 
theoretical approach. 
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Figure 14: The order parameter behavior within the first-order regime of the 
tricritical contact process (TCP). The solid line is obtained from steady state mea- 
surements at zero field for various system sizes. Within the simulations the control 
parameter p is slowly decreased until the absorbing state is reached. The circles 
correspond to data which are obtained from finite-field simulations and slowly in- 
creasing parameter p. Within the shadowed area the active and the absorbing 
phase coexist. The dashed line marks the transition point p c obtained from a sta- 
bility analysis of separated phases (see text). The inset shows the field dependence 
of the order parameter (for growing p). The three dashed lines correspond to three 
different field values (from h = 4 10 -4 (top) to h = 3 10 -5 ). As can be seen, a well 
defined upper limit of the supercooled low density phase can be obtained. 



6 FIRST-ORDER PHASE TRANSITION 

In this section we investigate the first-order regime of the modified contact process. 
A general phenomenon associated with first-order transitions is the presence of 
hysteresis by cycling across the transition. In equilibrium, the hysteresis is related 
to the effects of supercooling and superheating. Analogous effects occur in case of 
the modified contact process. First we investigate the effect of superheating, i.e., 
we consider the order parameter within the active phase (p a > 0) while approaching 
the transition point (on heating in equilibrium). Decreasing the parameter p for 
fixed q > q t the order parameter jumps at a certain point p from a finite value 
to zero. This is shown in Fig. El f° r Q = 0.97. Close to p = 0.221 the order 
parameter behaves discontinuously. This transition point p corresponds to the 
limit of superheating. Note that in contrast to a continuous phase transition the 
transition point p exhibits no systematic system size dependence. 

The similar phenomenon of supercooling the zero-order parameter phase is usu- 
ally not accessible for absorbing phase transitions. Owing to the lack of fluctuations 
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Figure 15: Snapshots of the tricritical contact process within the regime of first- 
order phase transitions (q = 0.97 and L = 256). The left figure shows a typical low 
density configuration (p = 0.223, supercooled phase). The seeds are triggered by 
an external field h = 10~ 4 . The right figure displays a high density configuration 
(p = 0.221, superheated phase). Here, order parameter fluctuations lead to seeds 
of various sizes. In both cases the seeds are subcritical, i.e., they disappear after a 
certain lifetime. 

for /? a = the system can never escape the absorbing state by a variation of p and q, 
respectively. To avoid that the system is trapped within the absorbing state an 
external field h is applied. In that way it is possible to perform steady state mea- 
surements of the order parameter on cooling, i.e., within in the low density phase 
(p a <C 1) for increasing p. The resulting curves are shown in the inset of Fig. 1141 At 
a certain value p u (h) the order parameter jumps from a low density value to a high 
density value. Note that the low density value of p a tends to zero for vanishing 
field. Furthermore, a well defined transition value p u (h) exists for h — > 0. The 
obtained value p n ~ 0.223 corresponds to the limit of supercooling. 

In that way we have obtained from steady state measurements a small but 
finite hysteresis, i.e., the two phases coexist between p < p < p u . Within the 
active phase (p a rs 0.8) the system is stable against small fluctuations until p 
is reached from above. On the other hand the absorbing phase is stable against 
external fluctuations, triggered by the conjugated field, until p u is approached from 
below. Snapshots of the system within the supercooled and superheated state are 
shown in Fig. 1151 

Finally we address the question of the critical value p c (q) of the first-order phase 
transitions. In equilibrium, the transition point is related to a thermodynamical 
potential such as the free energy. At the critical temperature the free energy of 
both phases are equal. Unfortunately, this definition can not be applied to the 
considered non-equilibrium phase transition. An alternative way of defining the 
first-order transition point is based on the behavior of moving interfaces which 
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Figure 16: Snapshots of the tricritical contact process close to the first-order 
transition for q = 0.97 and L = 128 (periodic boundary conditions are applied). 
Starting from a stripe of particles the dynamics is attracted either by the empty 
lattice (left, p = 0.22, the eventual occurring empty lattice is not shown) or by 
a steady state of a homogeneous non-zero particle density (right, p = 0.24). The 
time t is measured as the number of lattice updates. 

separate both phases. In case of phase equilibrium the interface velocity is zero 
whereas it is non-zero if one phase is favored by the dynamics. 

According to that picture we have investigated the phase propagation within 
the first-order regime. Initially the system contains a stripe (width L/2) of occupied 
particles. All lattice site outside the stripe remains empty. Depending on p and q, 
the system reaches after a transient either the absorbing phase or a steady state 
of a homogeneous non-zero particle density. Snapshots are shown in Fig. 1161 For 
q = 0.97 and L = 256 we have performed more than 50 runs for each value of p. The 
dynamics are attracted by the empty lattice in all runs for p < 0.2215. On the other 
hand the active phase is always approached for p > 0.2219. For p = 0.2217 both 
phases appear with probability of roughly 1/2. Thus an interval exists where both 
phases are favored by the dynamics. But we observe that this interval decreases 
with increasing system size L, i.e., a well defined transition point p c (q) exists in the 
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thermodynamic limit. For q = 0.97 the value p c = 0.2217 is obtained (see Fig. ll4|) 
yielding an asymmetric hysteresis between 0.221 < p < 0.223. 

The above analysis of the first-order transition is performed for a fixed value 
of q. Varying q the first-order transition line p c (q) as well as the borderlines of 
supercooling p u (q) and superheating p (q) can be determined. In order to limit the 
numerical effort we focus to the determination of p . Since the hysteresis is quite 
narrow (Ap/p c ~ 0.009) p u presents a sufficient approximation of the transition 
line. The corresponding values are plotted in Fig.0 Remarkably, the first-order 
line ends for q — > 1 at a finite p value in contrast to the mean field phase diagram 
(see Fig.Q}. 

7 CONCLUSION 

In summary, we have considered a modification of the well established contact pro- 
cess in order to study the process of tricritical directed percolation. Taking pair 
reactions into account the modified process exhibits a non-trivial phase diagram 
containing a tricritical point. The tricritical point separates a line of second-order 
phase transitions from a line of first-order phase transitions. The transition along 
the second-order line belong to the universality class of directed percolation. Per- 
forming a simple scaling analysis the tricritical point is determined with high ac- 
curacy. This allows a detailed analysis of the tricritical scaling behavior within the 
steady state. In particular, we have determined the tricritical exponents, universal 
scaling functions as well as universal amplitude ratios. The obtained values of the 
critical exponents as well as the universal amplitude ratios are listed in Tabled 
Additionally we have investigated the first-order regime. A hysteresis is found 
from steady state measurements. Owing to effects of metastability supercooling 
and superheating phenomena are observed. An analysis of the dynamical scaling 
properties of tricritical directed percolation will be published elsewhere. 

I would like to thank H.-K. Janssen, A. Hucht, G. Odor, and P. Grassberger 
for fruitful discussions. 

References 

[1] H. Hinrichsen, Adv. Phys. 49, 815 (2000). 
[2] G. Odor, Rev. Mod. Phys. 76, 663 (2004). 



[3] S. Liibeck, Int. J. Mod. Phys. B 18, 3977 (2004). 



Journal of Statistical Physics 123, 193 (2006) 



29 



[4] J. Marro and R. Dickman, Nonequilibrium phase transitions in lattice models 
(Cambridge University Press, Cambridge, 1999). 

[5] H.-K. Janssen, J. Phys.: Cond. Mat. 17, S1973 (2005). 

[6] H.-K. Janssen, Z. Phys. B 42, 151 (1981). 

[7] P. Grassberger, Z. Phys. B 47, 365 (1982). 

[8] S. Liibeck and R.D. Willmann, Nucl. Phys. B 718, 341 (2005). 

[9] H.-K. Janssen, Phys. Rev. E 55, 6253 (1997). 

[10] I. Jensen, Phys. Rev. Lett. 77, 4988 (1996). 

[11] A.G. Moreira, Phys. Rev. E 54, 3090 (1996). 

[12] R. Cafiero, A. Gabrielli, and M. A. Muhoz, Phys. Rev. E 57, 5060 (1998). 

[13] J. Hooyberghs, F. Igloi, and C. Vanderzande, Phys. Rev. E 69, 066140 (2004). 

[14] T. Vojta and M. Dickison, Phys. Rev. E 72, 036126 (2005). 

[15] E. Domany and W. Kinzel, Phys. Rev. Lett. 53, 311 (1984). 

[16] J.W. Essam, J. Phys. A 22, 4927 (1989). 

[17] S.S. Manna, J. Phys. A 24, L363 (1991). 

[18] M. Rossi, R. Pastor- Satorr as, and A. Vespignani, Phys. Rev. Lett. 85, 1803 
(2000). 

[19] S. Liibeck and PC. Heger, Phys. Rev. Lett. 90, 230601 (2003). 

[20] D. Zhong and D. benAvraham, Phys. Lett. A 209, 333 (1995). 

[21] J.L. Cardy and U. C. Tauber, Phys. Rev. Lett. 13, 4780 (1996). 

[22] T. Ohtsuki and T. Keyes, Phys. Rev. A 35, 2697 (1987). 

[23] T. Ohtsuki and T. Keyes, Phys. Rev. A 36, 4434 (1987). 

[24] H.-K. Janssen, M. Midler, and O. Stenull, Phys. Rev. E 70, 026114 (2004). 

[25] K.E. Bassler and D. A. Browne, Phys. Rev. Lett. 77, 4094 (1996). 

[26] F. Bagnoli, N. Boccara, and R. Rechtman, Phys. Rev. E 63, 046116 (2001). 

[27] A.P.F. Atman, R. Dickman, and J. G. Moreira, Phys. Rev. E 67, 016107 
(2003). 



Journal of Statistical Physics 123, 193 (2006) 30 

[28] T.E. Harris, Ann. Prob. 2, 969 (1974). 

[29] G. Odor, private communication , (2005). 

[30] R. Dickman and T.Tome, Phys. Rev. A 44, 4833 (1991). 

[31] R.M. Ziffand B.J. Brosilow, Phys. Rev. A 46, 4630 (1992). 

[32] B.J. Brosilow and R. M. Ziff, Phys. Rev. A 46, 4534 (1992). 

[33] S. Liibeck and R.D. Willmann, J. Phys. A 35, 10205 (2002). 

[34] H. Mori and K.J. McNeil, Prog. Theor. Phys. 57, 770 (1977). 

[35] S. Liibeck and R.D. Willmann, J. Stat. Phys. 115, 1231 (2004). 

[36] J.L. Cardy and R.L. Sugar, J. Phys. A 13, L423 (1980). 

[37] S.P Obukhov, Physica A 101, 145 (1980). 

[38] T. M. Liggett, Interacting particle systems (Springer, New York, 1985). 

[39] R. Dickman, Phys. Rev. E 60, 2441 (1999). 

[40] E. Luijten, H.W.J. Blote, and K. Binder, Phys. Rev. Lett. 79, 561 (1997). 

[41] S. Liibeck and P. C. Heger, Phys. Rev. E 68, 056102 (2003). 

[42] K. Binder and D. W. Heermann, Monte Carlo Simulation in Statistical Physics 
(Springer, Berlin, 1997). 

[43] S. Liibeck and H.-K. Janssen, Phys. Rev. E 72, 016119 (2005). 

[44] P. Grassberger, J. Phys. A 22, 3673 (1989). 

[45] J.C. Lee, J. Phys. A 24, L377 (1991). 

[46] C.A. Voigt and R.M. Ziff, Phys. Rev. E 56, R6241 (1997). 

[47] P. Grassberger and Y.-C. Zhang, Physica A 224, 169 (1996). 

[48] I. Jensen, Phys. Rev. A 45, R563 (1992). 



